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NOMENCLATURE 


symbol used to represent a gravitational term, 

A, = sin i cos 6 cos e - cos i sin 0 
1 a a a 

symbol used to represent a gravitational term, A„ = sin <t> cos 0 

d a a, 

symbol used to represent a gravitational term, 

A_ = cos i cos <t> cos 0 + sin i sin 0 

5 a a a 

lateral cyclic pitch measured from hub plane and in hub-body system, 
rad (or deg) 

components of total linear acceleration at the aircraft center of 
gravity expressed in the hub-body system, ft/sec 2 (m/sec 2 ) 

tip loss factor 

longitudinal cyclic pitch measured from hub plane and in hub-body 
system, rad (or deg) 

blade chord, ft (m) 

blade section drag coefficient 

blade section lift coefficient 

flapping damper coefficient, lb-ft-sec/rad (N-m-sec/rad) 
lag damper coefficient, lb-ft-sec/rad (N-m-sec/rad) 
center of gravity 

incremental drag force on blade element, lb (N) 
incremental lift force on blade element, lb (N) 
flapping hinge offset, ft (m) 
lag hinge offset, ft (m) 

distance between the flapping and lag hinges, ft (m) 
gravitational acceleration, ft/sec 2 (m/sec 2 ) 

tilt angle of the rotor shaft in the rotorcraft x a -z a plane, positive 
tilt forward, rad (deg) 
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2 

mass moment of inertia of the rotor blade about the lag hinge, slug-ft^ 
(kg-m 2 ) 

mass moment of inertia of the rotor blade about the flapping hinge, 
slug-ft 2 (kg-m 2 ) 

mass moment of inertia of the rotor blade segment between the flapping 
and lag hinges about the inner hinge (flapping or lag hinge), slug-ft 2 
(kg-m 2 ) 

pitch-flap coupling ratio 
pitch-lag coupling ratio 

flapping spring constant, lb-ft/rad (N-m/rad) 
lag spring constant, lb-ft/rad (N-m/rad) 

components of the offset between the aircraft center of gravity and the 
hub center in hub-body system (t x positive when the hub center is 
forward of aircraft c.g.; ly positive when the hub center is to the 
right of aircraft c.g.; *, z positive when the hub center is below 
aircraft c.g.), ft (m) 

mass of blade per unit of radius, slug/ft (kg/m) 

mass of blade from lag hinge to blade tip, slug (kg) 

mass of blade from flapping hinge to blade tip, slug (kg) 

local Mach number at blade element (M r = U/a, with a being the speed 
of sound) 

mass moment of rotor blade about lag hinge in a horizontal position, 
slug-ft (kg-m) 

mass moment of rotor blade about flapping hinge in a horizontal posi- 
tion, slug-ft (kg-m) 

mass moment of rotor blade segment between the flapping and lag hinge 
about the inner hinge (flapping or lag hinge depending on a hinge 
sequence) in a horizontal position, slug-ft (kg-m) 

components of aircraft angular velocity in hub-body system, rad /sec (or 
deg/sec) 

components of aircraft angular acceleration in hub-body system, 
rad/sec 2 (or deg/sec 2 ) 



generalized force associated with lead-lag for L-F-P hinge sequence, 
lb-ft/rad (N-m/rad) 

generalized force associated with flapping for L-F-P hinge sequence, 
lb-ft/rad (N-m/rad) 

generalized force associated with lead-lag for F-P-L hinge sequence, 
lb-ft/rad (N-m/rad) 

generalized force associated with flapping for F-P-L hinge sequence, 
lb-ft/rad (N-m/rad) 

generalized force associated with lead-lag for F-L-P hinge sequence, 
lb-ft/rad (N-m/rad) 

generalized force associated with flapping for F-L-P hinge sequence, 
lb-ft/rad (N-m/rad) 

aerodynamic force contribution to the generalized force associated with 
lead-lag and flapping, respectively (different expressions for the 
three hinge sequences), lb-ft/rad 

gravitational force contribution to the generalized force associated 
with lead-lag and flapping, respectively (different expressions for the 
three hinge sequences), lb-ft/rad 

contribution of dampers to the generalized force associated with lead- 
lag and flapping, respectively, lb-ft/rad 

contribution of hub springs to th^ generalized force associated with 
lead-lag and flapping, respectively, lb-ft/rad 

radial position of the blade element measured from the lag hinge, ft 
(m) 

• 

radial position of the blade element measured from the flapping hinge, 
ft (m) 

radius of blade root pocket cutout (i.e., radius at which lifting 
surface of blade begins), ft (m) 

rotor blade radius measured from the center of rotation, ft (m) 
total kinetic energy of a blade, lb-ft (N-m) 

resultant relative air velocity perpendicular to blade-span axis at 
blade element, ft/sec (m/sec), U = (u^ + Up) V 


u T component at the blade element of resultant relative air velocity 

perpendicular to the blade-span axis and on the plane parallel to 
x 1 -y 1 plane, ft/sec (m/sec) 

Up component at blade element of resultant relative air velocity perpen- 

dicular to both the blade-span axis and u^, ft/sec (m/sec) 

Ug component at blade element of resultant relative air velocity along the 

blade span, ft/sec (m/sec) 

V p total inertial velocity of the blade element, ft/sec (m/sec) 

6W fc total virtual work done on the blade by an external force, lb-ft (N-m) 

x , y , z hub-body system 

x',y',z' inertia system 

x a’ y a ,z a body-axis system 

i f y ,Z components of velocity at the aircraft center of gravity in hub-body 

system, ft/sec (m/sec) 

x,y,z components of velocity at blade element due to rotation, flapping, 

pitching, and lagging motion (i.e., isolated rotor) in hub-body system, 
ft/sec (m/sec) 

Ax,Ay,Az component of velocity at blade element due to aircraft motion in hub- 
body system, ft/sec (m/sec) 

X A ,Y A ,Z ft components of the total* aerodynamic force on blade element in hub-body 
system, lb (N) 

X ,Y r ,Z_ components of the gravitational force on blade element in hub-body 
system, lb (N) 

a r blade element local angle of attack, rad (or deg) 

<t> r tan" 1 (u p /u T ), rad (or deg) 

0 blade pitch angle measured from hub plane, 

0 = 0 o - A 1s cos * - B 1g sin + 0 fc r - - K 2 c, rad (or deg) 

0 Q blade-root collective pitch measured from hub plane, rad (or deg) 

0. blade twist at blade element, rad (or deg) 

t,r 

8 blade flapping angle measured from hub plane, rad (or deg) 


viii 



first and second derivative of 8 with respect to time 

blade lagging angle (measured positive when opposite to the direction 
of rotor rotation), rad (or deg) 

first and second derivative of c with respect to time 

Euler angles of the body-axis system with respect to the inertia system 
using the sequence of rotations as indicated 

mass density of air, slug/ft^ (kg/m^ 
rotor system angular velocity, rad/sec 

azimuth angle measured from -x axis in the sense of rotor rotation, 
rad (or deg) 



SUMMARY 


This report documents a derivation of coupled flap-lag equations of motion for 
a rigid articulated rotor with hinge springs and viscous dampers. Three different 
flapping-lag-pitch hinge sequences are considered and the Lagrange method is used in 
deriving the equations. The effects of the complete six degrees-of- freedom aircraft 
motions are included and all the inertia dynamic terms are retained; no small-angle 
assumptions are used in the development. Comparisons of the results of this work 
with those available in the literature are made. Sources of terms missing in previ- 
ous analyses, especially those of the inertia dynamics, are identified. 


INTRODUCTION 


Most helicopters in production today utilize articulated main rotors. A recent 
projection (ref. 1) indicates that, for the transport class of helicopters, future 
development may still rely heavily on the conventional articulated rotor design, 
although a trend seems to be set to use hingeless or bearingless rotors in the 
design of smaller combat helicopters to enhance agility and maintainability. 

An articulated rotor uses flapping, lead-lag, and feathering hinges (or elasto- 
meric bearings) in an appropriate sequence. The flapping hinge is employed to 
relieve the blade-root bending moment arising from the lift and centrifugal forces, 
the lag hinge to alleviate the in-plane moment at the blade root due to the drag 
force and the Coriolis force which stems from the flapping motion, and the feather- 
ing (or pitch) bearing to change the angle of attack of the blade. Figure 1 shows 
three widely used hinge sequences, i.e., flap-lag-pitch (F-L-P), flap-pitch-lag 
(F-P-L), and lag-flap-pitch (L-F-P) hinge arrangements. Use of a specific hinge 
sequence is often dictated by the past design experience and preference of the 
helicopter manufacturers. For example, Sikorsky designs tend to use an L-F-P hinge 
sequence with colocated hinges (refs. 2 and 3) as schematically shown in figure 2. 
The other two hinge sequences, the F-L-P and F-P-L, are more commonly used by other 
manufacturers. 

The flap-lag equations of motion and the dynamic characteristics of the articu- 
lated rotor blades have been under intensive investigations for a long time 
(refs. 4-11), largely by rotor dynamicists, but it is only relatively recently that 
the importance of the hinge sequence has been recognized (refs. 8 and 9). The rotor 
dynamicists tend to place their main emphasis on the stability characteristics of 
the isolated rotor system and, as a consequence, will tend not to include the 
effects of aircraft (or shaft) motions (refs. 4-9). Therefore, the flap-lag 
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HINGES: 1. FLAP 

2. LAG 

3. FEATHERING (PITCH) 



ROTOR ROTOR ROTOR 

SHAFT SHAFT SHAFT 

HINGE SEQUENCES: 


(F-L-P) (F-P-L) (L-F-P) 

Figure 1.- Three commonly used hinge sequences (ref. 2). 



ROTOR 

SHAFT 


Figure 2.- Elastomeric-bearing rotor hub system (ref. 2). 

equations of motion developed with this emphasis are not readily usable by the 
flight dynamicists, who are interested primarily in the stability and control char- 
acteristics of the entire vehicle and, as such, the effects of aircraft motions must 
be included in the development of flap-lag equations of motion. Some effects of 
aircraft motions were included, notably by the flight dynamicists (refs. 3, 10, 
and 11), but the effects are for specialized hinge arrangements. Extensions of 
those flap-lag equations of motion to include other hinge sequences and/or to 
include the complete effects of aircraft motion for a wider range of applications 
are not readily apparent, and therefore require a new development. 

In the past, most rotor models developed expressly for flight dynamics and 
control applications, especially those for use in real-time simulation, have consid- 
ered only the flapping motion (refs. 12-16). Recent studies related to high- 
bandwidth flight control systems for rotorcraft (refs. 17-21) have shown not only 
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that the flapping and other higher-order dynamics need to be considered 
(refs. 17-19), but the lead-lag dynamics must also be included in the rotor math 
model (refs. 20 and 21). These initial studies, coupled with the dramatically 
increased computational power that has become available in the past few years, serve 
as strong motivation to develop a generic rotor model for an articulated rotor 
system that includes various hinge arrangements for the rotor blades for a wider 
range of flight dynamics applications. It also permits a detailed investigation 
into the effects of hinge sequence. 

The work reported here represents the first of a series of efforts designed to 
achieve the above objective. A detailed derivation is provided herein in the devel- 
opment of flap-lag equations of motion for the three commonly used hinge sequences, 
i.e., L-F-P, F-P-L, and F-L-P for articulated rigid blade with a complete account 
for the effects of aircraft motion. 


COORDINATE SYSTEMS 


In this section, the position and velocity of an element of a rotor blade are 
calculated. To express these vector quantities, the following coordinate systems 
are used. Different coordinate systems are used for the three hinge sequences 
considered. 


L-F-P Hinge Sequence 

Figure 3 shows the coordinate systems that are used for the L-F-P hinge 
sequence. The blade is assumed to be rigid. The pitch bearing is outboard of the 
flapping hinge, and the flapping hinge is outboard of the lag hinge. 


SHAFT 



Figure 3.- Coordinate systems used for the L-F-P hinge sequence. 


The flapping hinge is assumed to be perpendicular to the rotor shaft. The lag 
hinge is assumed to be parallel to the shaft. The feathering axis is assumed to 
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coincide with the blade span axis, and the center of gravity and aerodynamic pres- 
sure center of the blade sections are assumed to be on the blade span axis. With 
these assumptions in mind, the following coordinate systems are used: 

x 1 ,y 1 ,z 1 - With its origin on the flapping hinge F; z 1 parallel to the shaft and 
pointing downward; y^ parallel to the blade segment between the two 
hinges (coinciding with the feathering axis when flap angle is equal to 
zero) and pointing outward; z y^, and x^ forming an orthogonal set. 

x 2 ,y 2 » z 2 ~ With origin on the lag hinge L; z 2 parallel to z^; y 2 parallel to the 
blade segment between the shaft and L and pointing outward; z 2 , y 2 , and 
x 2 forming an orthogonal set. 

x^.y^.z^ - With origin at the hub center, otherwise the same as x 2 ,y 2 ,z 2 system. 

This coordinate system will be referred to as the rotating shaft system. 

x,y,z - This axis system will be referred to as hub-body system with origin also 
at the hub center; z parallel to Zoj x pointing forward to the nose of 
the aircraft; y toward the right. This system is related to the rotat- 
ing shaft system through the azimuth angle ip as depicted in figure 4. 



Figure 4.- Relationship of the rotating-shaft system and the hub-body system. 


It is important to note that the sign convention regarding the flapping and lagging 
angles as defined in this report is: positive for flapping up and for lagging 

back. With this sign convention, the component of the position vector r with 
respect to the rotating-shaft and the hub-body systems are, respectively, 
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and 


r 



-(Ae + r cos 8)sin x. 

P 

+ (Ae + r cos b) cos c> 

P 

-r sin B 

B J 


r ■n 

x 

< y> 

z 


sin ip -cos ip 

cos <p sin ip 

0 0 



( 1 ) 


= < 


-(Ae + r 
-(Ae + r 


S 


B 


cos B)sin c sin 
cos B)sin c cos 


>P - [e^ + (Ae + r g 
ip + [e^ + (Ae + r g 
-r sin B 


cos B)cos c]cos >p 
cos 8)cos c]sin ip > 

j 


( 2 ) 


The hub-body system is now related to the aircraft body-axis system (x a ,y a ,z a ) and 
an inertia system (x',y' t z'). Figure 5 depicts the relationships among the three 
coordinate systems. Note that we consider a shaft tilt angle i, and a general 
offset of the hub center of the rotor system from the center of gravity of the 
aircraft (r H ). The center of gravity of the aircraft is translating (v Q ), acceler- 
ating linearly (a Q ), and rotating (u>,<I>) with respect to the inertia system. To 
expedite the development that follows, these vector quantities are now expressed in 
terms of the hub-body system, i.e., 

^ J 

r = (x,y,z) 


r u = U ,i ,l ) T 
H x’ y’ z 

v = (x ,y ,z ) 
o o* J o f o 


U) 


(p,q,r) 


T 


( 3 ) 


Note that when expressed in the hub-body system (also in the aircraft body-axis 

system), a = (d/dt)(v ) = (3/3t)(v ) + a> x v , and 
o o o o 
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(d/dt)u> = (3/3t)u> + <■> x at = (3/3t)u>, owing to the rotation of the system with 
respect to the inertia frame x',y',z'. Thus, 


a 


o 


<vv a z ,T 5 ( wV T * 


0 

r 

-q 


-r 

0 

P 


q 

-p 

0 


U,y,zr 


• t • • • \T 
U> = (p,q,r) 



Figure 5.- Relationships among inertia, body -axis, and hub-body systems. 


It should be noted that, due to the shaft tilt angle i, the components of the 
vector quantities in equation (3) differ from those which are expressed in the 
aircraft body-axis system by a transformation T^, 


“ cos 

i 

0 

sin 

i 

0 


1 

0 


-sin 

i 

0 

cos 

i. 


(4a) 


For example, if v q = (x ,y a ,z ) T and u> = (P a >q a » r a ) are, respectively, the 
linear velocity an§ angufar velocity of the aircraft c.g. when expressed in the 

aircraft body-axis system, then in the hub-body system, 

/A A A \ T m / A A A \T , / _ _ _ \ T T(n “ — 




a- z a> 


o’-'o 7 O l - _ _ 

sideward tilt with an angle 

is to be replaced by 

transformations . 


and (p,q,r) 


- *q *r 
a’ la 7 a 


If the shaft has a 


j, in addition to the forward tilt angle i, then 


T i 


Tj^ = T j T i’ whictl i- s a result of the sequence of rotational 
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The aircraft body-axis system, x a ,y a ,z a , is related to the inertia system, 
x',y',z', through the usual sequence of rotation of the Euler's angles, ^ • 

The components of a vector expressed in the two frames are related by the following 
transformation 




< 


> 


a 

v -/ 


cos 0 cos i|j 
a a 


sin 0 sin <t> cos 

a a a 

- cos 4> sin i|> 
a a 

sin 0 cos A cos ib 
a a a 

+ sin $ sin i|> 
a a 


cos sin A 
a a 


sin 9 sin $ sin i|i 
a a a 

+ cos $ cos ib 

a d 

sin 0 cos * sin i|> 
a a a 

- sin 0 cos ♦ 
a a 


_ 


r 

-sin 0„ 
a 


x' 

sin $ cos 0 


y’ , 

a a 

< 


cos 0 cos <b 
a a 


z' 

m 


V J 


(4b) 


From figure 5, it is clear that the sign convention for the components 
t x ,t y ,t z in equation (3) is as follows: t x is positive when the hub center H is 
forward of aircraft c.g. (along x-axis); ly is positive when H is to the right 
of c.g. (along y-axis); and fc z is positive when H is below aircraft c.g. 

(along z-axis). 


With the above considerations we are now in a position to determine the total 
velocity of the blade element. 


Total velocity of blade element - With reference to figure 5, the position 
vector with respect to the aircraft c.g. is: 


r r = r H + r 

T T 

Expressed in the hub-body system, r r = (t x + x, ! + y, i- z + z) , where (x,y,z) 

are given in equation (2). 

The total inertial velocity of the blade element V r is (fig. 5): 


V r = at <r r> ■ dt <'i> * 3t < r r> 


v o + it (r r ) + “ x r r 


(5) 


Note that (3/3t)(r r ) = (3/3t)(r), because (3/3t)(r^) = 0. When expressed in the 
hub-body system, (3/3t)(r) can be obtained by a simple differentiation of equa- 
tion (2) with respect to time to yield 
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t 


e a sin 4> + r e sin e cos(4> - c) + (Ae + r cos B)(a - c)sin(4> - c) 

€ p p 

e 8 cos 4> - r„B sin e sin(4> - c) + (Ae + r cos 6)(8 - c)cos(4> - c) 


-re cos e 

P 


(6) 


The other terms on the right-hand side of equation (5) are due to aircraft motion 


v + w x 
o 




-q 



po + qA • rB "1 

= < y n - pa + ^ > 


+ pB - qc 

v- ° -/ 


where from equation (2), 

A = a - r„ sin 6 
z 6 

B = (a + e_ sin 41) + (Ae + r cos e)sin(i|> 

y s d 

C = (a - e cos 41) - (Ae + r cos e)cos(4> 

X £ d 



(7a) 


(7b) 


(7c) 


In the subsequent development, the velocity components of the blade element 
(eq. (7b)), which are due to aircraft motions, will be denoted by (Ax,Ay,Az) . 
Equations (5), (6), and (7) thus completely define the total velocity components, in 
the hub-body system, of the blade element for the L-F-P hinge sequence. The coordi- 
nate systems used for the F-P-L hinge sequence are described next. 
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F-P-L Hinge Sequence 


We now consider the second hinge sequence, F-P-L. With the pitch bearing in 
board of the lag hinge, the blade pitch motion couples directly to the lead-lag 
motion of the blade, because the feathering axis no longer passes through the blade- 
span axis. A sequence of coordinate systems, which differs from that used for the 
L-F-P hinge sequence previously discussed, is now set up as shown in figure 6 . 

These coordinate systems have been previously used in reference 10 for a comprehen- 
sive flap-lag rotor analysis. That analysis was specifically for tandem rotor 
configurations. In order to apply the results of analysis to other rotorcraft 
configurations, such as side-by-side dual-rotor or tilt-rotor configurations, a more 
general analysis must be undertaken. 



ROTOR SHAFT 


Figure 6 .- Coordinate systems used for the F-P-L hinge sequence. 

In figure 6 , the x,y,z system is the hub-body system; the Xc,y^,Z(j system 
is the rotating-shaft system (shown as x^y^Zg in the L-F-P case); and the 
x 4 >Y 4» z 4 system is identical to the x^,y^,z^ system except that its origin is now 
shifted from the hub center H to point F on the flapping hinge line, which is 
assumed to be in a direction perpendicular to the shaft. The x^y^z^ system has 
its origin also located at F and is obtained by rotating about x^ (thus x^ = X 4 ) 
a flapping angle 8 so that y^ aligns with FL. The ^ 2 ^ 2' z 2 s y stem is 
obtained from x^y^z^ system by shifting its origin from F to L on the lag 
hinge. Finally, the x«,y«,z« system is obtained from the ^ 2 ^ 2 ,z 2 s Y stem b Y a 
rotation 0 about the y 2 axis. Now with respect to the x^,y 1 ,z 1 system, the 
blade element dm has the components: 
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( 8 ) 



With respect to the other coordinate systems, the components are obtained using the 
following transformations: 


and 



(9) 


( 10 ) 


( 11 ) 


(12) 


(13) 


Combining equations (8) through (13), the components of the position vector r from 
the hub center H to the blade element dr (with mass dm as shown in fig. 6), 
when expressed in the hub-body system, are obtained: 
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r = 


n 




~-(e + Ae cos 8 + r cos 8 cos c + r„ sin 8 sin c sin e)cos i|i 

8 c C 


- r^ cos 0 sin c sin 


= J (e 0 + Ae cos 6 + r cos 8 cos C + r sin 6 sin c sin 0)sin it» 
' p 5 v 

- r cos 9 sin c cos 41 
C 

-Ae sin 8 - r^(cos z. sin 8 - cos 6 sin x, sin 0 ) 


> ( 14 ) 


Total velocity of the blade element - As in the case of the L-F-P hinge 
sequence, the total velocity of the blade element dr is the sum of the relative 
velocity of the blade element with respect to the hub center and the velocity caused 
by the aircraft motion as given in equation (5). The relative velocity with respect 
to the hub center, when expressed in the hub-body system, is obtained as before, 
that is, by directly taking the derivative of equation (14) with respect to time, 
yielding 


St < r > ■< 


fa 


V- J 


(15) 
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where 


x 


fl[sin 4»( e^ + sin 8 sin 

- r^ cos i|> cos 0 sin c] + 

- r cos 8 sin z sin 0) + 
+ cr^(cos i|> cos 8 sin z - 


Z sin 0 + Ae cos 6 + cos 8 cos z) 

8 cos ♦( r* c sin 8 cos z + Ae sin 8 

0r^ sin c(sin ♦ sin 0 - cos sin 8 cos 0) 

sin <|» cos 0 cos z - cos i|> sin 8 cos z sin 0) 


y = fl[cos 4»(e + r sin 8 sin z sin 0 + Ae cos 8 + r 

p c c 

+ r sin 4» cos 0 sin c] + 8 sin *(-r c sin 8 cos z 

+ r cos 8 sin z sin 0) + 0r sin c(cos 4> sin 0 + 

C t 

+ cr^(-sin ip cos 8 sin z - cos ♦ cos z cos 0 + sin i|i sin 0 cos z sin 0) 

z = -8(r cos z cos 6 + Ae cos 8 + r sin 8 sin z sin 0) 

+ 0r^(cos 6 sin z cos 0) + cr^(sin z sin 8 + cos 6 cos z sin 0) 

When expressed in the hub-body system, the components of the velocity of the blade 
element due to the aircraft motion is given in equation (7), in which the components 
of the position vector, (x,y,z) T , are now given by equation (14) rather than equa- 
tion (2). Denote these velocity components by (Ax,Ay,Az) ; then 


cos 8 cos c) 

- Ae sin 6 

sin <|> sin 8 cos 0) 


(16) 


Ax 


V 

o 


0) X 


r 


r 


Ay 


> 


Az 


(17a) 
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where 


Ax = x + q{l - Ae sin 8 - r (cos c sin 8 - cos 8 sin c sin 0)} 

0 2 C 

- rU y + [e g + Ae cos 8 + r^cos 8 cos c + sin 8 sin c sin 8)]sin 

- cos 0 sin c cos 

Av = V + r{«, - [e + Ae cos 8 + r (cos 8 cos c + sin 8 sin c sin 6)]cos i|» 

J J o x 8 C 

- r cos 0 sin c sin i|>} - p{i - Ae sin 8 - r (cos c sin 8 

G z c 


(17b) 


- cos 8 sin x, sin 0)} 


Az = z - qU - [e + Ae cos 8 + r (cos 8 cos c + sin 6 sin c sin 9)]cos ij> 
ox6 C 

- r cos 0 sin x. sin >!>} + pU + [e + Ae cos 8 + r (cos 6 cos c 

C JT P 

+ sin 8 sin x. sin 0)]sin i> - r cos 0 sin c cos il>) 

and the components of the total velocity of the blade element are the sum of equa- 
tions (16) and (17), i.e., 


X + Ax 


V r = < y + Ay 


Z + AZ 


( 18 ) 


This completes the determination of the position and velocity components of the 
blade element in the hub-body system for the F-P-L hinge sequence. 


F-L-P Hinge Sequence 


The position and velocity of the blade element for the F-L-P hinge sequence can 
be obtained directly from those for the F-P-L hinge sequence by omitting the trans- 
formation caused by blade pitch rotation. Setting 0=0 in equation (14) yields 
the position components in the hub-body system: 


r 



/• 

-(e„ + Ae cos 
8 

< (e g + Ae cos 


8 + r^ cos 8 cos 

8 + r cos 8 cos 

-(Ae + r cos 
S 


c)cos ♦ 
c)sin t|> 
c)sin 8 


r 

C 


sin c 
sin c 


sin 

cos 



(19) 
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Likewise, the total velocity of the blade element is obtained by setting 
0=0=0 in equations (16) and (17b) to yield: 


V_ = < 


X + AX 
y + Ay 
z -i- Az 


(20a) 


where 


x = Q[sin ip(e + Ae cos 0 + r cos $ cos c) - r cos ip sin t] 

PS S 

+ 6 cos ip sin s( Ae + r^ cos c) + cr^cos ip cos 0 sin c - sin ip cos c) 


y = n[cos ip(e + Ae cos 0 + r cos 0 cos c) + r sin ip sin c] 

P S S 

- 0 sin ip sin 0(Ae + cos c) - cr^sin ip cos 0 sin c + cos ip cos c) 
z = -8 cos 8(Ae + cos c) + cr^ sin c sin 0 

and 


A x = x + q{» - (Ae + r cos c)sin 0} 

o z C 

- r{£ + [e. + (Ae + r cos c)cos 0]sin ip - r cos ip sin t) 

y p s s 


Ay = y + r{t - [e Q + (Ae + r cos c)cos 0]cos ip - r sin ip sin O 

0X65 C (20b) 

- p{& - (Ae + r cos c)sin 0} 

z X, 


Az = z - q{t - [e. + (Ae + r cos c)cos 0]cos ip - r sin ip sin e) 

O X p s £ 

+ p{t + [e + (Ae + r cos c)cos 0]sin ip - r cos ip sin ?} 

y 0 t t 

This concludes the derivation for the expressions of the position and velocity 
components of the blade element, in the hub-body system, for the F-L-P hinge 
sequence . 


The total velocity as given in equations (5), (18), and (20), respectively, for 
L-F-P, F-P-L, and F-L-P hinge sequences is required for developing the coupled flap- 
lag equations of motion using the Lagrange method. It is also required for the 
calculation of the aerodynamic forces and moments to be discussed later in this 
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report. The development of the flapping and lag equations is given in the next 
section. 


INERTIA DYNAMICS 


The Lagrangian procedure is applied to develop the inertia dynamics for the 
flapping and lag motion. This procedure begins with the determination of the 
kinetic energy of a single blade. This is obtained by integrating the kinetic 
energy of the blade element over the entire blade span. The development of flap/lag 
equations of motion will be given separately for each of the three hinge sequences. 


L-F-P Hinge Sequence 


The kinetic energy of the single rotor blade is determined by summing up that 
for each blade element, (1/2)V r • V r dm, over the entire blade span, i.e., 


T = 1/2 


r 

Jo 


V • V dm 
r r 


( 21 ) 


To help identify the contribution of the aircraft motion to the total kinetic 
energy, the dot product of the total velocity, V r • V r , will be expressed in the 
following manner: 

V p • V p = (x + Ax) 2 + (y + Ay) 2 + (z + Az) 2 


• 2 *2 * 2 \ "2 *2 • 2 * * • • •• 

=(x +y + z ) + Ax +Ay + Az + 2(Axx + Ayy + Azz) 


( 22 ) 


The first term of the right-hand side of equation (22) involves no effects of air- 
craft motion and it can be calculated directly from equation (6); the remaining 
terms involve the effects of aircraft motion and are determined using equations (6) 
and (7). The result is 

22 2*2 2 *2 
V * V = e_U + r 0 + (Ae + r Q cos 8 ; (Q - i) + 2e Q[r fl 0 sin 6 sin c 

r r g o p C p 

* 2 * 2*2 

+ (Ae + r cos 6 )(n - c ) cos c] + (xf + y„ + z*\ 

P Q O O ) 

+ 2 (qx - py )(t - r sin 6 ) + 2 r{y fU - e^ cos 4 >) 

- (Ae + r cos g)cos(ip - c)] - x [(j, + e sin i|>) 
p o y c 

+ (Ae + r cos 0 )sin(i|> - c)]> + 2z^{p[(t„ + e, sin 4 ») 
p o y c 

+ (Ae + r cos B)sin( 4 » - c ) ] - q[(t„ - e, cos 41 ) 

P X ^ 
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- (Ae + r cos 8)cos(H» - c)]} + 2{x [e Q sin if + r 6 sin 8 cos(H> - c) 

p o c 8 

+ (Ae + r cos B)(G - c)sin(t|> - c)] + y_[e o cos i|i 

P o ( 

- r a 8 sin 6 sin(t|> - c) + (Ae + r Q cos 8)(n - c)cos(i|> - c)3 

P P 

- 2 r0 cos 8} + (p + q )(* - r sin 8) 

Op Z p 

2 2 2 

+ (p + r )[(Ay + sin if) + (Ae + r g cos 6)sin(* - c ) ] 

2 2 2 

+ (r + q )[(«. - e, cos i|>) - (Ae + r cos 8)cos(i|> - c)] 

x c P 

- 2r(a - r sin 6){q[U„ + e r sin *) + (Ae + r 0 cos 8)sin(* - t)3 

z p y c p 

+ p((fc - e cos <|>) - (Ae + r cos 8)cos(i|» - c)]} 
x c p 

- 2pq[(fc + e sin ip) + (Ae + r cos 8)sin(ip - t)][(j, - e cos if) 

y c p x c 

- (Ae + r cos 8)cos(i|i - c)] 4- - r Q sin 8){q[e ft sin 4> 

P Z p C 

+ r 8 sin 6 cos(i|j - c) + (Ae ♦ r Q cos 6)(ft - c)sin(i|> - c)] 
p m 

- p[e ft cos \l> - r 8 sin 8 sin(\|> - c) 

C p 

+ (Ae + r cos 6) (0 - c)cos(* - c)]} - 2r 8 cos B(p[U„ + e r sin ip) 
p p y s 

+ (Ae + r g cos 6)sin(i|> - t) - qC ( «- x - e^ cos if) 

- (Ae + r cos 8)cos(ip - c ) ] } + 2r([(i - e„ cos *) 

- (Ae + r cos 8)cos(ip - c)][e Q cos if - r 8 sin 8 sin(<p - c) 

P Cm 

+ (Ae + r g cos 8) (a - c)cos(ip - c)] - [U y + sin i>) 

+ (Ae + r cos B)sin(<p - e)][e Q sin ip + r.B sin 8 cos(ip - c) 

P Cm 

+ (Ae + r g cos 6) ( Q - t)sin(ip - c ) ] > (23) 


The zeroth, first, and second mass moments of the blade segment outside the 
flapping hinge are denoted by 


I 


8 



dm , 



and 



dm 
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in which the limits of integration are from the flapping hinge to the tip of the 
blade. Consideration of the kinetic energy associated with the blade segment 
inboard of the flapping hinge introduces additional terms. Appendix A discusses the 
impact of this on the inertia dynamics for the three hinge sequences considered. 

The kinetic energy associated with a single rotor blade, T, is determined to be 

1 OOO o o o o 

T = 4 (8 I. + c ( Ae m 0 + 2AeM cos 8 + I fl cos 8) + 0 [(e^ + Aejm. 

2 2 

+ 2AeM cos 8 + I cos 8] + 20 e (Aem + M cos B)cos c 
p p C p 6 

oo 

- 2co[e (Aem + M cos B)cos c + (Ae m + 2AeM cos 8 + I. cos 8)] 

C 0 p p p p 

• 1 -2 *2 *2 • * 

+ 28Qe^Mg sin 6 sin c) + £ + Y 0 + z 0 ) m 6 + (qx Q - py o )U z m g - M g sin 8) 

+ r{y Q [U x - e^ cos 4»)m g - (Aem g + M g cos 8)cos(i|» - c)] 

* + e r sin *) m o + (Aem + M cos 8)sin(i|> - c)]> 

o y c p p p 

+ Z {pU + e sin 4>)m + (Aem + M cos B)sin(4» - c) - q[(l - e„ cos 4>)m„ 

vytppp x^p 

- (Aem Q + M cos 8)cos(<|> - c)]} + x re m a Q sin + M 0 8 sin 8 cos(i|> - c) 

+ (Aem g + M g cos 8)(0 - c)sin(4> - c)] + y o t e c m 6 ° cos * 

- M B sin 8 sin(i|» - c) + (Aem + M cos B)(Q - c)cos(<|> - c)] - z M 0 B cos 6 

P P p Op 

♦ \ (p 2 + q 2 )(^">g - 2a z Mg sin 8 + I g sin 2 8) 

+ \ (p 2 + r 2 )[U y + e c sin 4») 2 m e 

+ 2U y + e^ sin <i) (Aem g + M g cos B)sin(i|> - c) 

2 2 2 

+ (Ae m + 2AeM cos 8 + I. cos 8)sin ( tp - c)] 

P P P 

+ 5 (q 2 + r 2 )[U x - e c 003 <P) 2m 6 

- 2(l x - e^ cos 4» ) ( Aem g + M g cos 8)cos(i|> - c) 

2 2 2 

+ (Ae m g +• 2AeM g cos 6 + I g cos 8)cos (<l» ■* c)] 


- pq [ ( «. + e sin il>)U v - e cos 4»)m 

y q * <• p 

- («, + e sin i|>)(Aem + M cos 8)cos(i|» - c) 

y c p p 

+ U x - e^ cos il»)(Aem g + M g cos 8)sin(i|» - t) 

2 2 

- (Ae m g + 2AeM g cos 8 + I g cos 6)sin(i|> - c)cos(i(i - c)] 
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- pr(U - e cos 4>)U m - M sin 6) - [i (Aein + M n cos 0) 

X Q Zpp Zpp 

- AeM sin 0 - I sin 0 cos 0]cos(ip - c)} 

0 P 

- qr{(i y + sin <C)U z m g - M g sin 0) 

+ [i (Aem + M cos 0) - AeM sin 0 - I sin 0 cos 0]sin(4» - c)} 

Zpp p P 

+ (l m - M sin 0){q[e Q sin 4> + Ae(Q - e)sin(4i - c)] 

Zpp £ 

- pte^Q cos 4> + Ae(Q - c)cos(4> - ?)]} 

+ Ot z M 6 - I g sin 0){q[0 sin 0 cos(4> - c) + (0 - c)cos 0 sin(4> - c)] 

+ p[0 sin 0 sin(4> - c) - (0 - c)cos 0 cos(i|i - c)]} 

- 0 cos 0 {p [ ( 8, + e sin 4i)M + (AeM + I cos 0)sin(i|» - c) 

y c o 8 p 

- q[U - e cos )M - (AeM + I cos 0)cos(i|» - ?)]} 

X £ p P p 

+ r{(l - e cos i|i)[m e a cos 4 ) - M o 0 sin 0 sin(i|» - 5 ) 

+ (Aem g + M g cos 0 )(fl - c)cos(4> - ?)] 

- ( 8 . + e sin 4 »)[m e 0 sin 4 » + M 0 sin 0 cos(4> - t) 

y c p c p 

+ (Aem + M cos 0 )(a - c)sin(\|» - c)] - (Aem + M cos 0)e a cos c 

P P P P C 

- (Ae 2 m g + 2AeMg cos 0 + I g cos 2 0 )(a - £)} (24) 

It can be seen from equation (24) that the aircraft motion greatly complicates the 
expression for the kinetic energy of the rotor blade. It also complicates, as can 
be expected, the equations of motion of the blade. 

The coupled flap-lag equations of motion are now derived using the Lagrange 
procedure (refs. 22 and 23). They are calculated using 



(25a) 



(25b) 


In equations (25a) and (25b), and Q g are the generalized forces (moments in the 
present case). They consist of the effects due to gravity, aerodynamic forces, and 


18 



the forces due to lag dampers and hinge springs. These will be discussed later in 
this report. 

Coupled flap-lag equations of motion - With some algebraic manipulations of 
equations (24) and (25a) the lag equation of motion becomes: 

"2 2 
c(Ae m. + 2AeM 0 cos 8 + I Q cos 8) 

DP P 

2 2 

= 0(Ae nig + 2AeMg cos 8 + Ig cos 8) + (Aenig + Mg cos 8)ta x sin(ip - c) 

+ Sy cos(ip - t) + e[(6 - r)cos 5 - (0 - r) sin z - pq cos(2ip - c) 

- p sin ip cos(\p - c) + q cos ip sin( 4> - c)] - Ae[r + pq cos 2(ip - z) 

+ ^ (p 2 - q 2 )sin 2(4* - t)l + & x [(r + P^) 003 ^ - c) - (q^ + r 2 )sin(ip - t)] 

+ t y [(pq - r)sin(<p - c) - (p + r )cos(ip - c)] 

- t [ (p - qr)cos(* - c) - (q + pr)sin(ip - c)]} 

z 

+ (AeM + I cos 6) {sin 8[2B(r - 0 + c) + (p - qr)cos(ip - z) 

6 D 

12 2 

- (q + pr)sin(4> - c)] - cos 8tr + pq cos 2(<p - z) + 2 (P * Q ) sin 2 (<P ~ 5) 

- 26[p cos(ip - c) - q sin(<p - c ) ) > } + Q r (26) 

Similarly, from equations (24) and (25b), the flapping equation of motion is 
obtained: 

- M a • 2 

8 = 7^ {sin B(a sin(<p - c) - a cos(ip - z) - e[6 sin t + (0 - 2r0)cos c] 

i 0 y x 

- Ae[(Q - C) 2 - 2r(Q - z)] 

* 2 2 

+ [t - e cos - Ae cos(ip - c)][(r + pq)sin(<p - z) + (q + r )cos(ip - c)] 

X c 

• 2 2 

+ [£y + e^ sin tp + Ae sin(ip - c)][(r - pq)cos(ip - z) - (p + r )sin(ip - c)] 

- - qr)sin(4> - z) + (q + pr)cos(ip - z)]} 
z 

+ cos eta + 2p[oe cos \p + Ae(Q - c)cos(\p - c)] 

z r z 

- 2q[Qe j , sin ip + Ae(Q - c)sin(tp - c)l 

- (q - pr) [1 - e cos ip - Ae cos(ip - c)] 

X (» 

2 2 

+ (p + qr ) [ t + e_ sin ip + Ae sin(\p - c)] - (p + q H }} 
y s ^ 
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+ sin e{[(p - qr)sin(\|i - c) + (q + pr)cos(i|» - c)]sin 6 

2 *2 

- { [p sin(4> - c) + q cos(i|> - t)] + [r - (Q - c)] }cos 8} 

2 2 * 

+ cos B { ( p + q )sin 8 + {2(fl - c)[p cos(i|> - t) - q sin(\|» - c)] 

+ (p + qr)sin(4» - c) + (q - pr)cos(i|> - c))cos b} + t — (27) 

ft 


In equations (26) and (27), a x , ay, and a z are, respectively, the components of the 
total linear acceleration at the aircraft c.g. when projected along the x,y,z axes 
of the hub-body system. It is important to recognize that these components are 
different from x Q , y Q , and z Q , which are obtained by differentiating x Q , y , and 
z q with respect to time, without considering the fact that the hub-body system is 
rotating about the inertia system, x',y',z'. In fact, from equation (3), 

a x = x 0 * p y 0 + qi 0 (28a) 

a y = y 0 + rk O " P^O (28b) 

a 2 = z 0 + py 0 - qi 0 (28c) 

This distinction becomes particularly important when combining the flap-lag equa- 
tions with the rigid-body aircraft equations of motion to calculate the dynamics and 
responses of the complete rotorcraft. 

Equations (26) and (27) are the coupled flap-lag equations of motion for the 
L-F-P hinge sequence. It should be noted that in deriving these equations, the 
inertial effects of the blade segment inboard of the flapping hinge were 
neglected. When these effects are included, the inertia dynamic terms in the lag 
equation (26) must be expanded to include additional terms given by equation (A2) in 
appendix A. However, the flapping equation (27) remains unchanged, because they 
have no direct effect on the flapping motion. 


F-P-L Hinge Sequence 

The coupled flap-lag equations of motion for the F-P-L hinge sequence are 
developed in the same manner as for the L-F-P hinge sequence discussed previously. 
Using equations (16), (17b), (21), (22), (25a), and (25b), the coupled flap-lag 
equations of motion are obtained: 
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= 8 sin 6(1 + AeM cos c) ♦ fl cos © [ I cos 0 + M cos c(e + Ae cos 8)] 


+ 2001^ cos 0 sin z 


- 200 cos 0 cos c[I^(sin 0 cos z - cos 0 sin z sin 0) + AeM^ sin 0] 

m 2 

+ 2Q0 sin cl^(sin 0 cos c - cos 0 sin c sin 0) + 0 1^ sin c cos £ 

• P 2 

- (j sin c(I cos c cos 0 + AeM ) 


+ a 


[2 sin 2c ^ sin2 6 " cos2 8 sir)2 0 ) + cos 2c sin 0 cos 8 sin ej 


+ (e.M + AeM cos 0)(sin 0 cos c sin 0 - cos 6 sin c)V 
B € t J 

+ M fa (cos c sin 4* cos 0 + cos c cos * sin 0 sin 0 - sin c cos 4> cos 0) 

C x 

+ a (cos c cos * cos 0 - cos c sin 4> sin 0 sin 0 + sin c sin ip cos 6) 

y 

- a (cos c sin 0 cos 0 + sin c sin 0)] 

z 

+ 20 1 {sin c cos c(q sin 4> cos 0 - p cos ip cos 0 - r sin 8) 

+ sin^ cfsin 0(q sin 4> sin 8 - p cos 4> sin 8 + r cos 8) 

- cos 0(q cos tp + p sin 4>) ] > 

+ 28[(I t cos z cos 0 + AeM^sin c(q cos * + p sin ip) 

- (1^ cos z + AeM^Jcos z cos 0(q sin 4> cos 8 - p cos 4> cos 8 - r sin 8) 

- I cos 0 sin 0 cos z sin c(q sin 4> sin 8 - p cos 4> sin B + r cos 8)] 

+ 20{f(e M + AeM cos 8) + I (cos 8 cos z + sin 0 sin z sin 0)] 

0 c c t 

x [sin 0 cos c(q sin ip cos 8 - p cos ip cos 8 - r sin 0) 

+ sin c(q sin 4* sin 8 - p cos ip sin 6 + r cos 8)] 

- I cos 0 sin c[(sin 0 cos z cos 6 + sin z sin 8)(q cos 4> + P sin ip) 

+ r cos z cos 0]} + {M [-ft + (e + Ae cos 6)cos ip] 

£ X D 

+ I^[sin <p cos 0 sin z + cos ip(cos 8 cos z + sin 8 sin c sin 0)]} 

OO • • 

x {sin 0 cos cl(q + r )cos ip sin 0 + (pq + r)sin ip sin 8 - (q - pr)cos 8] 

- sin c[(q + r )cos 4> cos 8 + (pq + r)sin 4> cos 6 + (q - pr)sin 8] 
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2 o 

+ cos £ cos 0[(q^ + r )sin ip - (pq + r)cos ip]} 

+ {M^y + < e g + Ae cos 6)sin ip] 

+ I^[-cos 0 sin c cos ip + sin ip(cos 6 cos x, + sin 0 sin £ sin 0)]} 

2 2 * 

x {sin 0 cos c[(pq - r)cos ip sin 0 + (p^ + r )sin ip sin 0 - (p + qr)cos 0] 

• 2 2 

- sin c [ ( pq - r)cos ip cos 0 + (p + r )sin ip cos 0 + (p + qr)sin 0] 

2 2 

+ cos 0 cos c[(pq - r)sin ip - (p + r )cos ip]} 

+ {M (-Z + Ae sin 0) + I (sin 0 cos x, - cos 0 sin x, sin 0)} 
t z s 

x {sin 0 cos t[-(q + pr)cos ip sin 0 - (p - qr)sin ip sin 0 

o o 

- (p^ + q^)cos 0] - sin c[-(q + pr)cos ip cos 0 — ( p - qr)sin ip cos 0 

2 2 § * 

+ (p + q ;sin B] + cos 0 cos c[-(q + pr)sin \|> + (p - qr)cos il>]} + Q (29) 


2 2 2 2 
B[Ae m + 2 AeM cos c + I (cos c + sin c sin 0)] 

= c sin 0(1^ + AeM^ cos c) + (0 cos 0 sin x. + 2c0 cos 0 cos c - 0 sin 0 sin c) 

x (I cos c + AeM ) + (200 sin 0 sin c - 20c cos 0 cos c - 0 cos 0 sin c) 

£ £ 

x [I (sin 0 cos 0 sin x, - sin 0 cos £) - AeM sin 0] 

2 

+ 2&i sin c( 1^ cos £ cos 0 + AeM^) 

- 2001^ sin 0 cos 0 sin^ £ - i AeM^ sin 0 sin c 


4* 


2 2 2 2 
sin 20[Ae + 2AeM^ cos c + I^(cos x, - sin 0 sin c)] 


+ 2 ^ sin cos 20 sin 0 + AeM^ sin 0 sin x, cos 20 

+ e.M (cos 0 sin x, sin 0 - sin 0 cos c) - Ae e.ra sin 0 
p C p c 


) 


- a cos ip[m Ae sin 0 - M (sin 0 cos 0 sin c - sin 0 cos £)] 

x c C 

- sin ip[-m^ Ae sin 0 + M^(sin 0 cos 0 sin x. - sin 0 cos £)] 

+ a [-ra Ae cos 0 + M (sin 0 sin 0 sin x. + cos 0 cos c)] 

- 2c{-(q sin ip sin 0 - p cos ip sin 0 + r cos 0)1^ sin 0 cos 0 sin x, cos £ 
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+ (q sin ip cos 8 - p cos ip cos 6)cos 0 cos c(M^ Ae - 1^ cos c) 

+ r sin 8 cos 0 cos c(M^ Ae + cos c) + (q cos ip + p sin ip) 

x [M^ Ae sin c(1 - 2 cos 2 8) + 2 Ae sin 6 cos 6 cos c sin 0 

2 

+ 1^ sin x. cos c cos 0]> 

2 . 2 

- 20{(q sin ip sin 6 - p cos ip sin 8 + r cos 8)1^ sin 0 sin c 
+ (q sin ip cos 8 - p cos i> cos 8)sin 0 sin c(-M^ Ae + 1^ cos c) 

- r sin 0 sin c sin 8(M^ Ae + 1^ cos c) 

+ (q cos ip + p sin ip)sin c cos 0(2M^ Ae sin 8 cos 8 - sin c sin 0)} 

2 

-2s(q cos ip + p sin ip)[2M^ Ae sin 0 sin c cos 8 

- 2 Ae sin 8 cos 8(M t cos c + in, Ae)] 

-2Q{sin 0 sin c(-p cos ip sin 6 + q sin ip sin 8 + r cos 8) 

x [M e„ + I (sin 8 sin z sin 0 + cos 8 cos c)] 

- sin 0 sin z cos 8 Ae(-p cos ip sin 8 + q sin ip sin 8 - r cos 6) 

- (p cos ip cos 6 - q sin ip cos 8 + r sin cos c(e g + Ae cos 6) 

+ Ae sin 8 sin 0 sin z + cos t(cos 8 cos z + sin 8 sin z sin 0)] 

+ (p cos ip cos 8 - q sin ip cos 8 - r sin B ) [m^ Ae(e g + Ae cos 8) 

+ M Ae cos 8 cos t] - cos 0 sin c(p sin ip + q cos ip) 

c 

x [-M^ Ae cos 8 + I^sin 0 sin z sin 8 + cos z cos 8)]) 

2 2 

+ [(q + pr)cos ip sin 6 + (p - qr)sin ip sin 8 - (p + q )cos 6] 

x (m Ae(-fc + Ae sin 6) - M Ae(cos 8 sin z sin 0 - cos z sin 8)) 

£ z c 

2 2 

+ [(q + pr)cos ip sin 8 + (p - qr)sin ip sin 8 + (p + q )cos 8] 

x (M cos c(-i + Ae sin 8) + I cos c(cos C sin 8 - sin z sin 0 cos 8)> 

C z C 

2 2 

- [(q + pr)cos ip cos 6 + (p - qr)sin ip cos 8 - (p + q )sin 8] 

x {M sin 0 sin c(-l + Ae sin 6) 

C Z 
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+ I sin z sin 0(sin 8 cos c - cos 0 sin z sin 0)> 

- [(r - qp)cos ip sin 0 - (p + r )sin 4> sin 0 - (p + qr)cos 0] 

x (M^ Ae[cos ip cos 0 sin z - sin 4>(cos 0 cos c + sin 0 sin 0 sin c)] 

- m Ae(e fl + Ae cos 0)sin 4»> 

C P 

o o 

- [(f - pq)cos 4> sin 0 - (p + r )sin 4* sin 0 + (p + qr)cos 0] 

x {-a m Ae - M cos c[a + (e + Ae cos 0)sin i(>] 

y z z y b 

+ 1^ cos c[cos 4i cos 0 sin c - sin 4>(cos 0 cos z + sin 0 sin 0 sin t)]} 

- [(r - pq)cos i|) cos 0 - (p + r )sin * cos 8 - (p + qr)sin 0] 

x (M^ sin 0 sin + (e g + Ae cos 6)si n ♦! “ sin Z s ^ n 6 

x [cos 4> cos 0 sin z - sin 4»(cos 0 cos z + sin 0 sin 0 sin c)]> 

+ [(q 2 + r 2 )cos 4i sin 0 + (pq + r)sin 4> sin 0 + (q - pr)cos 0] 

x (m Ae[J,„ - (e 0 + Ae cos 0)cos 41] 

C x p 

- Ae[sin 41 cos 0 sin z + cos 4i(cos 0 cos z + sin 0 sin 0 sin c)]} 

+ [(q 2 + r 2 )cos 41 sin 0 + (pq + r)sin 4» sin 0 - (q - pr)cos 0] 

x (M cos t[»„ - (e a + Ae cos 0)cos 4»] - I r cos z 

C x p *-» 

x [sin 41 cos 0 sin z + cos 4)(cos 0 cos z + sin 0 sin 0 sin t)]} 

+ [(q 2 + r 2 )cos 41 cos 0 + (pq + r)sin 4» cos 0 + (q - pr)sin 0] 


x {-M sin 0 sin t[t - (e + Ae cos 0)cos 4>] + I sin 0 sin z 

t X p 4 

x [sin 41 cos 0 sin z + cos 4i(cos 0 cos z + sin 0 sin 0 sin z)D + Q g 


(30) 


In equations (29) and (30), m^, , and 1^ are the zeroth, first, and second 

mass moments of the blade segment outboard of the lag hinge, i.e., 

•R /*R 

1 . 1 

m 


p r 

Jo z Jo 


r dm , I 

z * z 


f 


r 2 dm 

z 


The limits of integration for these quantities are from the lag hinge to the tip of 
the blade. When the effect of the inertia dynamics of the blade segment inboard of 
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the lag hinge is considered, additional terms as given by equation (A4) must be 
added to the right-hand side of the flapping equation (30) as discussed in detail in 
appendix A. 


F-L-P Hinge Sequence 

Now with the pitch bearing outboard of the lag hinge, the motion of the blade 
pitch no longer couples to the flap and lag motion. The coupled flap-lag equations 

of motion for the F-L-P hinge sequence can therefore be obtained directly from 

equations (29) and (30) by setting e = e = e = 0. Alternatively, they can also be 
obtained using equations (20a), (20b), (21), (22), (25a), and (25b). The results 
are: 

I r C = 0[I_ cos 6 + M cos c(e 0 + Ae cos 8)] - 208 cos c sin e(I cos c + AeM ) 

C G C p C 5 

# P 0 f* i o 

-8 sin c(I t cos c + AeM^) + or I ^ 1^ sin^ 8 sin 2c 

- M (e + Ae cos 8)cos 8 sin c + M [a (cos c sin i]» - sin c cos ip cos 8) 

C p j c x 

+ a (cos c cos ip + sin c sin ip cos 6) - a sin c sin s] 

y 2 

+ 28(I C cos c + AeM^Hsin c(q cos ip + p sin p) 

- cos c(q sin ip cos 8 - p cos ip cos 6 - r sin 6)] 

+ 20{(M (e + Ae cos 8) + I cos 8 cos dsin c 
C p c 

x (q sin ip sin 8 - p cos $ sin 8 + r cos 6) 

- I sin c[(q cos ip + p sin dsin c sin 6 + r cos cl) 

+ (M [ - Jl + (e + Ae cos 8)cos <p] + I [sin ip sin c + cos ip cos 8 cos cl) 

C X fi £ 

2 2 

x {-sin ct(q + r )cos ip cos 8 + (pq + r)sin ip cos 8 + (q - pr)sin 8] 

p p 

+ cos c t ( q + r )sin ip - (pq + r)cos ip J } 

+ [ fty + (e g + Ae cos 8)sin ip ] - I^[cos ip sin c - sin ip cos 6 cos d) 

x {-sin t((pq - r)cos ip cos 8 + (p^ + r )sin <p cos 8 + (p + qr)sin 8] 

+ cos c[(pq - r)sin ip - (p^ + r )cos p]> 

+ (M (-ft + Ae sin 8) + I sin 8 cos d 

£ Z £ 
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X 


2 2 . 

{-sin c[-(q + pr)cos ip cos 8 - (p - qr)sin ip cos B + (p + q )sin B] 
+ cos c [ — ( q + pr)sin ip + (p - qr)cos </]> + Q 

**o 


(3D 


"2 2 
B[ Ae + 2 AeM^ cos x, + cos c] 


(28c cos c + 8 sin c)(I cos c + AeM )sin B + 26c sin c( I cos c + AeM ) 


4 5 ' 


o 2 

(Ae + 2 AeM^ cos x, + 1^ cos c)sin 26 


- e g sin B(M t cos x. + Aem^l - sin 6(M ? cos x, + Aera^)(a x cos ip - a y sin ip) 


+ cos 8 (M cos x, - Aem )a - 2c(cos c[cos B(q sin i|> - p cos ip) 
C Z 

x (M Ae - I cos O + r sin b(M Ae + I cos c)l 

C X. x, x. 

2 

+ sin £(q cos + p sin 4>)[M Ae(1 - 2 cos $) + I cos t]} 

+ 2B(q cos ip + p sin ip)(M cos c + m Ae)Ae sin 28 


+ 28 (cos c(p cos ip cos 6 - q sin ip cos B + r sin 6) 


x [M (e + Ae cos 6) + I cos 6 cos c] + sin c(p sin ip + q cos ip) 

C B C 

x [-M^ Ae cos B + 4 cos 5 cos 

- (p cos i> cos B - q sin tp cos 8 - r sin b) 

x [m Ae(e. + Ae cos 8) + M Ae cos 6 cos cl) 

2 2 

+ [(q + pr)cos 41 sin B + (p - qr)sin ip sin B - (p + q )cos B]Ae 

x [m (-a, + Ae sin 8) + M cos x, sin 6] 

1 x, z X, 

2 2 

+ [ ( q + pr)cos ip sin B + (p - qr)sin ip sin B + (p + q )cos B]cos c 

x [M (-«, + Ae sin 8) + I cos x, sin B] 

C Z Q 

o o • 

_ [(f - qp)cos ip sin B - (p + r )sin ip sin B - (p + qr)cos B]Ae 

x [M (cos ip sin c - sin ip cos 8 cos c) - m (e + Ae cos B)sin ip] 

2 2 * 

- [(r - qp)cos ip sin 6 - (p + r )sin ip sin B ♦ (p + qr)cos 6] 

x {-a m Ae - M, cos + (e a + Ae cos B)sin ip] 

y c z y p 

+ I cos c[cos ^ sin c - sin ip cos 6 cos c]) 
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0 5 » a 

+ t(q ♦ r )cos ♦ sin 8 + (pq ♦ r)sin ♦ sin 6 + (q - pr)cos 8]Ae 

x {m [a - (e + Ae cos 8)cos *] - M [sin $ sin c + cos i|» cos 8 cos *;]} 

C x p t 

2 2 * * 

+ [(q^ + r )cos i|> sin 8 + (pq + r)sin i|> sin 8 - (q - pr)cos 8]cos e 

x (M (a - (e + Ae cos 8) cos ♦ ] - I [sin i|» sin c + cos 4* cos 6 cos c]} 

C x 0 C 

+ Q 0 (32) 

S 3 

Equation (31) provides the complete inertia dynamics for the lead-lag equation 
of motion. However, as with the F-P-L hinge sequence, additional terms, as given 
by equation (A4), must be included in the right-hand side of the flapping equa- 
tion (32), when the effect of the inertia dynamics associated with the blade segment 
inboard of the lag hinge is considered. 


COMPARISON OF THE RESULTS WITH PREVIOUS WORK 


The inertia dynamics in the coupled flap-lag equations of motion for the three 
hinge sequences, as shown in equations (26) through (32), are now compared with 
available results in the literature. In general, the previous results have tended 
to be only for one specific hinge sequence or with limited considerations for the 
effects of aircraft motions. 


L-F-P Hinge Sequence 

A complete derivation of the inertia dynamics for the coupled flap-lag equa- 
tions of motion, including some effects of aircraft motions, is given in refer- 
ence 11 for an L-F-P hinge sequence with co-located flap and lag hinges. Refer- 
ence 3 also provides an approximation for the same configuration. This specific hub 
configuration can be obtained by setting Ae = 0 in equations (26) and (27) to 
yield: 


M. 


t = 


(a x sin(4» - s) + a cos(i|> - c) 


I cos 6 l “x ' ~y 

P 


+ e^Ua - r)cos c - (Q - r ) c sin t - pq cos(2ip - x.) 

2 2 

- p sin * cos(<) - c) + q cos sin(i|> - c)] 

• 2 2 

+ * x t(r + qp)cos(t|> - c) - (q + r )sin(i|> - c)] 

2 2 

+ ft y [ ( pq - r)sin(il) - c) - <p * r )cos(4> - s)3 
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- I [(p - qr)cos(i|> - c) - (q + pr)sin(* - c)]} 

z 

+ tan g[26(r - Q + £) + (p - qr)cos(4> - c) - (q + pr)sin(i|> - c)] 

+ (Q - r) + 2B[p cos(i|> - c) - q sin( 4 i - c)] - pq cos 2(4/ - c) 

- \ (P 2 - q 2 )sin 2 ( 4 * - O + (33) 

I cos B 

P 

and 

- M a • 2 

6 = j {sin 6 (-a cos(ip - c) + a sin(* - e) - e [n sin t + (Q - 2rfl)cos cl 

1 0 x y * 

® • 22 

+ (i - e cos t|») [ (r + pq)sin(i|/ - c) + (q + r )cos(i|> - t)] 

2 2 

+ U + e sin i|>)[(r - pq)cos(4> - c) - (p + r )sin(4> - ?)] 

...J r ** 

- «- [ (p - qr)sin(4> - c) + (q + pr)cos(4» - c)]} 

Z 

+ cos eta + 2ae (p cos 4 * — q sin i|>) - (q - pr)(Jl - e cos 41 ) 

Z 4 x £ 

• 2 2 
+ (p + qr)(t^ + sin 4») - ( P + q H z >} 

+ sin e{ [ (p - qr)sin(4> - c) + (q + pr)cos(4» - c)]sin B 

2 *2 

- { [p sinU - c) + q cos(4» - c)] + [r - (a - c)] }cos b} 

2 2 * 

+ cos b{(p + q )sin 6 + {2(a - c)[p cos(4> - c) - q sin(4» - t)] 

\ 

+ (p + qr)sin(i|» - t) + (q - pr)cos(4i - t)}cos b} + j — (34) 

*4 n 


Note that the underlined quantities in equations (33) and (34) are not included in 
reference 11 because the effects of the displacement between the hub center and the 
aircraft c.g. have not been considered. These effects are considered in the present 
development, allowing applications to tandem or side-by-side dual-rotor configura- 
tions where the offset between the hub center and aircraft c.g. are generally large 
and therefore cannot be neglected. Note also that the sign convention used in this 
report for the lead-lag angle is positive for lagging (roughly, measured positive 
when opposite to the direction of rotor rotation). Several other reports (e.g., 
refs. 3 t 5, 9, and 11) use different sign convention. 
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F-P-L Hinge Sequence 


A comprehensive flap-lag rotor analysis was reported in reference 10 for the 
F-P-L hinge sequence for tandem-rotor applications. The present work expands the 
effects due to aircraft motion to account for lateral offset !y which is necessary 
for modeling tilt-rotor or side-by-side dual-rotor configurations. The inertia 
moment for the blade segment inboard of the lag hinge, which was neglected in refer- 
ence 10, is also included in the present work for completeness (appendix A). 


F-L-P Hinge Sequence 


This hinge sequence (as well as the L-F-P hinge sequence) was considered in 
reference 9 for centrally co-located flap and lag hinges; no hinge offsets were 
included. Nor were the effects of aircraft motions considered in that work. Refer- 
ence 5 provides a list of the inertia dynamic terms for this hinge sequence. Unfor- 
tunately, neither the effects of the variation in the angular velocity of the rotor 
system nor the effects of aircraft motion were included. As discussed earlier, the 
effects of aircraft motions, which need to be considered for flight dynamics and 
control applications, drastically complicates the coupled flap-lag equations of 
motion as evident from equations (3D and (32). When these effects are neglected, 
equations (31), (32), and (A4) in appendix A reduce to simply 


I c : Q[I cos $ + M (e„ + Ae cos B)cos c] - 2fl8 sin 8 cos c(I cos c + AeM ) 
£ £ £ p 5 C 


• 2 c 

6 sin c(I cos c + AeM ) + Q 4 
£ £ 


I sin 2 B sin 2x. - M (e + Ae cos 6)cos 6 sin + Q 

c c p J s 


(35) 


2 2 

[Ae m + 2 AeM cos c + I cos t + I. ]B 
t C C Ae 

= (2ac cos c + a sin c )(I_ cos c + AeM ) sin 6 + 2 bc sin c(I cos c + AeM ) 

£ C £ ** 

2 Tl 2 2 

- a ~ Ue m + 2 AeM cos c + I„ cos c)sin 2 b + e. sin b(M cos c + Aem ) 

t £ 0 £ £ 

. <I 4e oo 3 6 * e B M 4e )sin s] ♦ Q 


(36) 


When the effects of the variation in rotor rpm, as indicated by the underlined terms 
in equations (35) and (36), are neglected, these equations reduce to those shown in 
reference 5. 

We now turn to examine two special cases: (1) co-located hinges with another 

hinge sequence, F-L-P, and (2) elimination of the lead-lag degree of freedom. 

When the flapping and lag hinges are co-located, the inertia dynamics for the 
coupled flap-lag equations can also be readily obtained by simply setting Ae = 0 
in equations (31) and (32). It is important to note, however, that the resulting 
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expressions will be different from equations (33) and (34) for the L-F-P hinge 
sequence, because the sequence of rotation associated with flapping and lag angles 
is different, thereby resulting in a different definition associated with these two 
angles (similar to the situation when the sequence of rotation is changed, the 
definition of Euler angles must change accordingly). It is for this reason that 
equations (26) and (27) were used as the basis for equations (33) and (34). 

Finally, in the absence of the lead-lag degree of freedom, it can be seen from 
equations (26) through (32) that the inertia dynamics for all three hinge sequences 
collapse to the same identical flapping dynamics as shown in equation (37). 

I 8 = -fl^feM sin 8 + 4 I sin 26j - [(a cos <l> - a sin 4>)sin 6 - a cos 8]M 
B y B do / x y z B 

+ 2Q(p cos ip cos 8 - q sin i|» cos 8 + r sin 6)(M e„ + I cos 6) 

Bp B 

* * 2 2 
+ [(q + pr)cos i|> sin 8 + (p - qr)sin p sin 8 + (p + q ;cos 6] 

• 2 2 

x (-MJ, + I sin 6) + [(r - pq)cos * sin 8 - (p + r )sin il> sin 6 

B Z p 

+ (p + qr)cos 8][M (ft + e sin ip) + I sin <l> cos 8] 

B y B s 

2 2 • • 

+ [ ( q + r )cos ip sin 8 + (pq + r)sin sin 6 - (q - pr)cos 6] 

x ^ M e^x ~ e s cos ~ 003 * cos + Q 6 

Equation (37) represents an extension of reference 14 by including all the flapping 
effects due to aircraft motion. 


DETERMINATION OF THE GENERALIZED FORCES 


The generalized forces (refs. 22 and 23) associated with the flapping and lead- 
lag motions, which are symbolically shown on the right-hand side of equations (26) 
through (32), consist of contributions from the gravitational force, aerodynamic 
forces, and the forces due to lag dampers and hinge springs. Determination of those 
contributions requires first calculating the virtual work of an individual contrib- 
uting external force associated with the virtual flapping and lead-lag displace- 
ments. Let X^YjjZj be the components of the ith external force F^ acting on 
the blade element dr when expressed in the hub-body x,y,z frame of reference. 
Then the resulting elemental virtual work done by this external force due to the 
virtual flapping and lag displacements 68 and 6c is 

6VL = X^6x + lL6y + Z^6z (38a) 
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where 


V * H ** * H " <38b) 


and x,y ,2 are the coordinates of the blade element as given by equations (2), 
(14), and (19) for the three hinge sequences being considered. Now, summing up the 
elemental virtual work over the appropriate blade span results in the total virtual 
work 6W{. 
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i 3t J 


6C 


(39a) 


By equating the total virtual work 6W t to the virtual work of the general i 2 ed 
forces along the virtual displacement in the generali 2 ed coordinates, i.e., 

6W t = ^6 56 + the generali 2 ed forces associated with that external force Fj 

are determined: 



(39b) 


(39c) 


As an example, the generali 2 ed forces Q_ and Q„ in equations (31) and (32) 

c 3 *3 

for the F-L-P hinge sequence will be derived in detail below to illustrate the 
procedure . Let 


Q 


C 


3 


Q 


6 


3 


= Q . 

♦ Q A 

+ Q _ 

+ Q c 

c,G 

C,A 

C,D 

c,s 

= Q 6,G 

+ Q 6,A 

+ Q 6,D 

+ Q 6,S 


(40a) 

(40b) 


representing the contributions from the gravity force, aerodynamic forces, and the 
forces due to dampers and hinge springs. We will consider each of these forces 
separately . 
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Gravity Force Contribution, q and Q g q 

The gravity force acting on a blade element with mass dm can be expressed in 
the hub-body system of reference using a sequence of transformation, equations (4b) 
and (4a), from inertia to aircraft body to hub-body systems. The result is 


X Q = gdm A 1 
Y g = gdm A 2 
Z Q = gdm Ag 

where 


(41) 


A. = sin i cos <t> cos 0^ - cos i sin 8 

1 cl el 3. 

A~ = sin <t> cos 0 (42) 

c a cL 

A- = cos i cos <t> cos 0 + sin i sin 0 

3 a a a 

From equation (19), the virtual displacement, in the hub-body system, of the blade 
element outboard of the lag hinge is obtained. 

6x = (Ae sin B + r^ sin B cos c)cos ip68 + r^(cos B sin c cos ip - cos c sin ip)6c 

6y = -(Ae sin B + r^ sin B cos c)sin ipAB - r^(cos 0 sin c sin ip + cos t cos ip)6c 

6z = -(Ae + r cos c)cos 868 + r^ sin c sin 86c 

(43) 


The coordinates for a blade element between the two hinges, measured r g from the 
flapping hinge, are 


x = -(e + r cos 0)cos ip 

3 3 

y = (e + r cos B)sin ip (44) 

3 3 

z = -r sin B 
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and therefore, 


fix = r sin 6 cos i|)60 
3 

fiy = -r sin 0 sin ipfig 

M 

fiz = -r cos 060 

P 

which are independent of the virtual displacement in the lead-lag motion. 

Substituting equations (41) —(43) and (45) into equation (39), the desired 
generalized forces due to gravity are obtained. 


(45) 


Q ? G = -gM^[(A^ sin cos i|/)cos x. + (Ag sin <> - cos i|»)cos 0 sin e 


- A^ sin 0 sin c] 


(46) 


Qg G = gUem^ + cos c + M a e ^ A 1 cos ^ " A 2 sin 8 - A^ cos 6] (47) 

Equations (46) and (47) represent a generalization of the results given previously 
in references 5 and 10 by including in Q g both the effects of shaft tilt angle 

and the weight moment inboard of the lag hinge. Note that equations (46) and (47) 
can alternatively be obtained using appropriate coordinate transformations listed in 
appendix B, and by considering the blade weight moment about the lag and flapping 
hinges as were done in references 5 and 10. 


Aerodynamic Force Contribution, Q g A and Q ? A 

The elementary aerodynamic forces acting on the blade element, dL and dD, are 
resolved first in the plane perpendicular to the blade span in a manner as illus- 
trated in figures 7(a) and 7(b). 


where 


dD. = dD cos $ - dL sin 4 

1 r r 

dL. = dL cos $ + dD sin $ 

i r r 


dD 

dL 


1/2 pU 2 c.c dr 
a t 

1/2 P U 2 c c dr 
® c 

2 2 
+ Up 

tan -1 (u p /u T ) 


(48a) 

(48b) 

(49a) 

(49b) 

(490 

(49d) 
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Figure 7.- Elemental aerodynamic forces, (a) Force diagram in x.|-y.j plane; 
(b) resolution of elementary lift and drag forces. 


Note that the section lift and drag coefficients c^ and c^ in equations (49a) 
and (49b) are functions of the local angle of attack a r and Mach number M r 
(ref. 24). 1 The calculation of a r and the airspeed components u p and u^. will be 
discussed later. 

When expressed in the x^y^z-j frame of reference, the components of the 
aerodynamic force are 



(50) 


From equations (B7) and (B8), these components can be transformed into those in the 
hub-body system 



-dD^cos c sin i|> - sin c cos 8 cos ip) + dL^ sin 8 cos ip I 
-dD^cos c cos <|> + sin e cos 6 sin ip) - dL 1 sin 8 sin ip > (51) 

dD 1 sin c sin 6 - dL^ cos 8 I 


Using equations (51), (43), (39b), and (39c), the aerodynamic force contributions to 
the generalized forces in equations (40a) and (40b) are obtained. 


1 

'Extensions of the two-dimensional airfoil data may be made to include the 
spanwise flow effect (refs. 25 and 26), unsteady effects (refs. 26-29), and dynamic 
inflow (ref. 30). 
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£t 


pU c,c sin $ (Ae + r cos c)dr 
d r C c 


Ci 


pU c c cos $ (Ae + r cos c)dr 
x r c C 


r 


1 2 

x pU c .c cos $ r dr 

2d T r c C 


-r» 


pU c.c sin $ r dr,. 
8. r c C 


(52) 


(53) 


Note that in setting up the limits of integration in equations (52) and (53), the 
blade cutoff is assumed to extend beyond the blade segment outboard of the lag 
hinge. Note also that in calculating the total blade lift, a tip loss factor B is 
used, thereby effectively limiting the lift of the blade beyond BR. 


The calculation of the components of the velocity (with respect to the air 
mass), Up and u T , and the local angle of attack, a r , is now addressed. The velocity 
of the blade element due to blade and aircraft motions has been obtained previously 
as shown in equation (20). Additional contributions to the total local airspeed 
components in the hub-body system are due to rotor-induced velocity and wind and 
turbulence. The resultant airspeed components in the hub-body system are then 
transformed into x 1 ,y 1 ,z 1 frame of reference using equation (B9). Let the resul- 
tant local airspeed components be denoted by x^,y^,z.. Then, u p , u T , and u s (the 
spanwise component) are given by 


Up = (54a) 

u T = x 1 cos t + y. sin x, (54b) 

u s = -x y sin t + y. cos c (54c) 

The local angle of attack of the blade element consists of 4> r and the blade 
pitch setting, which includes, for generality, collective and cyclic pitch, blade 
twist, and pitch-flap and pitch-lag couplings, and K 2 , as shown in equation (55). 


“r = *r + 9 o " A 1s cos * " B 1s sin * + 0 t,r " K 1 0 " K 2 C (55) 

This concludes the formulation of quantities that are necessary to compute the 
aerodynamic contribution to the generalized forces. 
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Contributions from Spring Restraints and Dampers 

The generalized forces corresponding to the spring restraints and dampers can 
be obtained directly (ref. 23) from the potential energy for the hub springs, 

(1/2)K B 2 + (1/2)K c , and the dissipation functions for the dampers, 

•2 ***2 

( 1/2)C B + (1/2)C^c , respectively. For generality, a provision is made, permit- 
ting either the representation of the actual springs/viscous dampers or the simula- 
tion of the structural stiffness /damping. These generalized forces are: 


Vs ■ -V 

(56a) 

Vs * -v 

(56b) 

Vd ■ -v 

(56c) 

Vd - c b® 

(56d) 


This concludes the development of all the factors contributing to the general- 
ized forces in equations (40a) and (40b) for the F-L-P hinge sequence. The general- 
ized forces for the other two hinge sequences considered can be developed in a 
similar fashion. Results are shown in appendix C. 


CONCLUDING REMARKS 


Flap-lag equations of motion for a rigid articulated rotor blade have been 
developed for three widely used flap, lag, and pitch hinge sequences. The effects 
of the complete six degrees-of-freedom of aircraft motions were included and all the 
inertia dynamic terms were retained. Comparisons of the results of this develop- 
ment, which uses the Lagrangian method, with those available in the literature were 
made to identify, with special emphasis on the inertia dynamics, the source of terms 
missing in previous analyses. 

This effort represents the first step toward developing a generic blade-element 
rotor model suitable for full-flight envelope, real-time, flight dynamic simulation 
of various rotorcraft. Subsequent steps will involve the development of detailed 
equations for the calculation of the rotor forces and moments, and the development 
of appropriate approximate models for specific applications. To permit comparisons 
with flexible blade analysis planned for a later date, a provision is made in the 
present analysis to allow for representation of the actual hub springs and viscous 
dampers or simulation of the structural stiffness and damping in both the flapping 
and lead-lag degrees of freedom. In addition, kinematic feedback couplings in 
pitch-flap and pitch-lag have also been included. 


I 
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APPENDIX A 


INERTIA DYNAMICS OF THE BLADE SEGMENT BETWEEN THE FLAPPING AND LAG HINGES 


In the development of the coupled flap-lag equations of motion, the inertia 
dynamics associated with the blade segment between the flapping and lag hinges were 
neglected. For completeness, additional terms due to this blade segment are devel- 
oped here for the three hinge sequences considered. 


L-F-P Hinge Sequence 


For this hinge sequence, the portion of the rotor blade inboard of the flapping 
hinge involves only the lead-lag motion. As a result, the inertia dynamics associ- 
ated with this segment contribute only to the lead-lag equation of motion (26); 
however, they have no effects on the flapping equation of motion (27). 


The total velocity of a blade element between the flapping and lag hinges, a 
distance r^ from the lag hinge, can readily be obtained from equations (6) and (7) 
by setting 0=0=0, Tg=0, and by replacing Ae with r^. The results are: 


where 


V = 


r. 

x + Ax 


<y + Ay 

+ Az 


x = e 0 sin * + r (a 
y = e^Q cos if + r^(fl 


c)sin(4» - c) 
c)cos(4» - t) 


(Ala) 


(Alb) 


z = 0 


and 

Ax = x + qt - r[(* + e sin \|>) + r sin(4» - ?)] 

o z y ^ £ 

Ay = y n - pt_ + r[(i - e cos i|>) - r cos(4» - c)] 

O Z X Q 

Az = z + p[(fc + e sin i|>) + r sin(4> - c)] 
o y c c 

- q[(£ - e cos i>) - r cosU - c)] 

X Q £ 


(Ale) 
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Using equations (21) through (25a), the inertia dynamics associated with the blade 
segment can be obtained as shown in equation (A2). These additional terms are to be 
added to the right-hand side of equation (26) to provide the complete set of inertia 
dynamics for the single rotor blade. Note that equation (A2) can also be directly 
obtained from equation (33) by setting 6 = 8 = 0 and by replacing M g and I g 
with M Ae and I Ae , respectively. 


•I Ae c + I Ae Q 6 - r) - ^ (p 2 - q 2 )sin 2(4> - t) - pq cos 2(<|> - t)] 

• • 2 

+ M {a sin(<|> - c) + a cos(i|> - c) + e [(0 - r)cos 5 - (0 - r) sin c 

ag x y £ 


- pq cos(2i|> - c) 



sin i|> cos(i|i - 


C) 


2 

q 


cos ^ sin(i|i - t)] 


where 


2 2 

+ * x [(r + pq)cos(i|» - c) - (q + r )sin(ip - c)] 

+ a y [(pq - r)sin(ip - z) - (p + r )cos(i|> - c)] 

- I [ ( p - qr)cos(4> - t) - (q + pr)sin(4> - 5)]} 

z 



(A2) 


F-P-L and F-L-P Hinge Sequences 

The effect of the inertia dynamics of the blade segment inboard of the lag 
hinge on the flapping equation of motion is identical for the F-P-L and the F-L-P 
hinge sequences. This blade segment involves only flapping motion and therefore its 
inertia dynamics have no effects on the lead-lag equation of motion, equation (29) 
or (31). 

The total velocity of a blade element between the flapping and lag hinges, 
measured r g from the flapping hinge, can be obtained from equations (20a) 
and (20b) by setting t = c = 0 and by replacing Ae with r g to yield equa- 
tions (Ala), (A3a), and (A3b). 

x = (e + r cos 8)8 sin i|> + r 8 cos sin 8 
SB S 

y = (e + r cos 6)0 cos <t> - r 6 sin sin 6 (A3a) 

SB S 

z = -r 0 8 cos 6 

p 


1 
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(A3b) 


Ax = x + qU - r sin 8) - r[a + (e + r cos 8)sin i|»] 
o z 0 y 0 p 

Ay = y Q + r[l x - (e g + r g cos 8)cos ip] - p(* 2 - r g sin 8) 

Az = z - q[a - (e. + r cos 6)cos ip] + pU + (e + r cos 6)sin ip] 

o x p p y p p 

Using equations (21) through (24) and (25b), the inertia dynamics associated with 
the blade segment are determined as shown in equation (A4). The same expression can 
also be obtained from equation (37) by simply replacing M g with M Ae , and I„ 
with I Ag . Adding these additional terms to the right-hand side of equations (30) 
and (32) provides the complete set of inertia dynamic terms for these two hinge 
sequences. 

-W ' ” 2 ( e 6 M 4e sln 6 * 5 ! 4e sln 2S ) 

- [(a cos il» - a sin *)sin 8 - a cos 8]M 
x y z 

+ 2fl(p cos <p cos 8 - q sin * cos e + r sin 8)(e 6 M Ae + 1 Le cos 6) 

+ [(q + pr)cos i(i sin 8 + (p - qr)sin ip sin 8 

+ (p 2 + q 2 )cos e](-M Ae 4 z + J Ae sin B) 

+ [ ( f - pq)cos i|i sin 8 - (p + r )sin ip sin 8 + (p + qr)cos 8] 

x [M. (a + e. sin i|>) + I._ sin it cos 8] 

1 Ae y 6 Ae 

+ [(q 2 + r 2 )cos ill sin 8 + (pq + r)sin it sin 6 - (q - pr)cos 6] 

x [M Ae ( a x - e g cos il») - I Ae cos 41 cos 6] (A4) 
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APPENDIX B 


USEFUL TRANSFORMATION MATRICES 


A list of useful transformation matrices is given in this appendix to facili- 
tate the determination of the generalized forces in the Lagrange equations of 
motion, (25a) and (25b). Specific transformations of interest vary with the hinge 
sequences considered. 


L-F-P Hinge Sequence 

The frames of reference of particular interest for this hinge sequence are: 
x^y^z,; x 3>y3» z 3i and x,y,z. The components of an arbitrary vector a in these 
frames are related by 




% 


cos c -sin x, 0 


sin c cos c 0 


£ ^ 


(B1 ) 






' a yf = 




sin ip -cos ♦ 0 


cos ip sin ip 0 


0 


0 1 


x 3 




(B2) 


J V U 




< a 

A 


y ■ 


V J 


sin ip cos x, - cos ip sin c cos ip cos c + sin ip sin x, 0 


-sin ip sin x. - cos ip cos c -cos ip sin x, + sin ip cos x, 0 


■\ 


<V (B3) 


a 

z 

v. J 
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Since these transformations are orthogonal, the inverse transformations are obtained 
by a simple transposition of the transformation matrices. 


F-P-L Hinge Sequence 


The frames of interest are: x^y^z^; X3, 73,23; and x,y,z. However, these are 

defined differently from the L-F-P hinge sequence (fig. 6). The components of a 
vector a in these three frames are related by the following transformation. 


r a ^ 
x 3 




% 

V. 3 


cos 0 0 sin e 


-sin 0 0 cos 0 






v>l 


(B 4 ) 


r 

a 

x 






sin ip -cos ip cos 8 -cos ip sin 6 


cos ip sin ip cos 8 sin ip sin 6 


-sin 8 


cos 6 


a 




(B 5 ) 






fsin ip cos 0 


V J 


cos ip cos e 


cos sin 8 sin 0 - sin ip sin 8 sin e 


-cos ip cos 6 


sin ip cos 6 


-cos 8 sin 0' 


-sin 6 


sin ip sin 0 cos ip sin 0 cos 8 cos 0 

L - cos ip sin 6 cos 0 + sin ip sin 6 cos 0 


/a^ 


( b6 ) 


2 

J 


F-L-P Hinge Sequence 

The frames of interest are the same as those for F-P-L hinge sequence. How 
ever, with the pitch bearing now outboard of the lag hinge, the transformations 
reduce to: 
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(B7) 


(B8) 


(B9) 
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APPENDIX C 


GENERALIZED FORCES FOR L-F-P AND F-P-L HINGE SEQUENCES 


A detailed derivation for the generalized forces in the coupled flap-lag equa- 
tions of motion for the F-L-P hinge sequence has been given earlier in the report. 
This appendix provides a summary of the results for the other two hinge sequences, 
i.e., L-F-P and F-P-L, using the same procedure for the F-L-P case. 


L-F-P Hinge Sequence 

Similar to equations (40a) and (40b) for the case of F-L-P hinge sequence, the 
generalized forces and Q g ^ in equations (26) and (27) also consist of four 

contributing sources, i.e., gravity force, aerodynamic forces, and forces from lag 
dampers and hub springs. The last two sources are the same as equation (56) for the 
case of F-L-P. The first two sources differ markedly from equations (46), (47), 
(52), and (53) as shown below. 

Gravity force contribution - The virtual displacement is developed from equa- 
tion (2) to yield 


6x = r sin 8 cos(i|> - c)6B - (Ae + r„ cos 8)sin(4» - c)6c 

P P 

6y = -r sin 8 sin(i|> - c)66 - (Ae + r a cos e)cos(i|» - c)6c (Cl) 

P p 

6z = -r cos 668 

P 

for a blade element outboard of the flapping hinge, and 

6x = -r sin( 4* - c)6c 

dy = -r cos (4» - c)6c (C2) 

6z = 0 


for inboard of the flapping hinge. Substituting equations (41), (42), (Cl), 
and (C2) into (39), the gravity force contribution to the generalized forces 
and Q_ is obtained. 


1 


Q 8,G = gM 8 {[A 1 C0S( * - C) 


A 2 sin(4» - c)]sin 6 - A^ cos 8} 


(C3) 


Q c,G = + M g cos 6 + M Ae^ A l sin ^' 1 ’ - c) + a 2 cos (i|) - e)] 
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Contribution of the aerodynamic forces - With reference to figure Cl, the ele- 
mentary lift and drag, dL and dD, as given by equation C5), 

dD = (1/2)pU 2 c d c dr g 


dL = (1/2)pU c^c dr g 

2 2 2 
U 2 = u 2 + u 2 

♦ r = tan” 1 (u p /u T ) 


(C5) 


are first resolved in the plane perpendicular to the blade span, 

dD^ = dD cos 4> r - dL sin $ r 

dL, = dL cos + dD sin $ 

1 r r 

and then in the x 1 ,y 1 ,z 1 frame of reference as shown in equation (C7) 


(C6) 



-dD, 


-dL^ sin B 
-dL 1 cos B 


(C7) 



These aerodynamic force components are further resolved into the hub-body system to 
yield 
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(C8) 


l: 


Z A 

v. *V 


-dD.j sin(i|> - c) + dL 1 cos(tp - c)sin 0 
-dD^ cos(4> - c) - dL.j sin(i|> - c)sin 0 
-dL.j cos 0 


Using equations (C8), (Cl), (39b), and (39c), the aerodynamic force contributions to 
the generalized forces and Q are obtained. 


. pi 

Jr 2 


pU c d c sin (d r r g dr g 


n 


pU c.c cos $ r„ dr,, 
l r 0 0 


(C9) 


- P 1 

Jr 2 


pU c d c cos ♦ r (Ae + r g cos 0)dr g 


r» 


pU c c sin <|> (Ae + r cos 8)dr 
* r 0 B 


(CIO) 


We note in passing that the transformation of the local airspeed at the blade ele- 
ment to u T , u s , and u p , shown previously in equation (54) for the F-L-P case is now 
replaced by 



The computation for the local angle of attack remains the same as 
hinge sequence using equations (55) and (49d). 


(Cl 1 ) 


that for the F-L-P 


F-P-L Hinge Sequence 

The contributions due to the lag dampers and the hub springs are the same as 
equation (56), but the other two factors become somewhat more complicated than 
equations (46), (47), (52), and (53) as shown below. 



Gravity force contribution - The virtual displacement is obtained from equa- 
tion ( 14 ) for both outboard and inboard segments of the lag hinge. Proceeding in 
the same manner as before for the F-L-P hinge sequence results in equations (C12) 
and ( C 1 3 ) . 

Q. r = g[(Aem + M cos c + M A )sin 6 - M sin x. cos 8 sin 0] 
x (A 1 cos i|> - A 2 sin i>) 

-g[(Aem + M cos x. + M )cos 8 + M sin x. sin 8 sin 0]A~ (C12) 

t c Ae c j 

Q r = gM {(cos 8 sin c - sin 6 cos x, sin ©HA. cos - A ? sin *) 

- cos 0 cos c(A^ sin + A 2 cos i|>> + (sin 8 sin c + cos 8 cos c sin 0)Ag} 

( C 1 3 ) 


Aerodynamic force contributions - Using the same procedure as before for the 


F-L-P case, the aerodynamic force contributions to Qg^ and 


are found to be 


Q 


8, A 



- sin 


1 2 

x pU c c[cos * cos 0(Ae + r cos c) 
c. Xt r 

4 sin 0(r + Ae cos c)]dr 

r C C 

1 2 

2 pU c d c[cos $ r sin 0(r^ + Ae cos c) 

<t> cos 0(Ae ♦ r cos c)]dr 
r z c 


- fi 

•« J, 2 


pU c .c cos <|>„r 


dr 

r c C 


-jn 


pU c^c sin 4 r r dr^ 


(Cl 4 ) 


(C15) 


The transformation (54) remains valid for computing the u p , u^., and ug from 
the local airspeed at the blade element in frame of reference. Equa- 

tion (B6) in turn provides the necessary transformation of the local airspeed compo- 
nents from the hub-body system to x^,y^,z^ system. 


The computation of the local angle of attack requires a special treatment for 
the F-P-L hinge sequence, because in this case 4» r specifically contains the blade 
pitch 0 and its time derivative 0, as evident from equations (16) and (17). It 
should be noted also that the inertia dynamics for this hinge sequence require 
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values of 6 and 6 as input quantities. Let 0 v be the root collective pitch 

when the lag hinge is parallel to the z axis. Then e is given by (ref. 10) 


0 = e 0 * 6 vpo * A 1s 003 * ' B 1s sin * ’ K l e - K 2 c 

Thus, in general, 

0 = 0 O - (A 1s + B 1s Q)cos * - (B ls ' A 1s Q)sin ^ ’ K i 0 - 
Finally, the local angle of attack is: 

a = <b + 0 +0. 

r Y r vpo t,r 


(Cl 6 ) 


( C 1 7 ) 


( C 1 8 ) 
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